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Abstract 

We investigate the possibility that the dark matter consists of clusters of the heavy family quarks 
and leptons with zero Yukawa couplings to the lower families. Such a family is predicted by the 
approach unifying spin and charges as the fifth family. We make a rough estimation of properties 
of baryons of this new family members and study possible limitations on the family properties due 
to the direct experimental and the cosmological evidences. 



I. INTRODUCTION 



Although the origin of the dark matter is unknown, its gravitational interaction with the 
known matter and other cosmological observations require from the candidate for the dark 
matter constituent that: i. The scattering amplitude of a cluster of constituents with the 
ordinary matter and among the dark matter clusters themselves must be small enough, so 
that no effect of such scattering has been observed, except possibly in the DAMA/Nal pQ 
and not (yet?) in the CDMS and other experiments [2]. ii. Its density distribution (ob- 
viously different from the ordinary matter density distribution) causes that all the stars 
within a galaxy rotate approximately with the same velocity (suggesting that the density 
is approximately spherically symmetrically distributed, descending with the second power 
of the distance from the center, it is extended also far out of the galaxy, manifesting the 
gravitational lensing by galaxy clusters), iii. The dark matter constituents must be stable 
in comparison with the age of our universe, having obviously for many orders of magnitude 
different time scale for forming (if at all) solid matter than the ordinary matter, iv. The 
dark matter constituents had to be formed during the evolution of our universe so that they 
contribute today the main part of the matter ((5-7) times as much as the ordinary matter). 

There are several candidates for the massive dark matter constituents in the literature, 
known as WIMPs (weakly interacting massive particles), the references can be found in [HE]. 
In this paper we discuss the possibility that the dark matter constituents are clusters of a 
stable (from the point of view of the age of the universe) family of quarks and leptons. Such 
a family is predicted by the approach unifying spin and charges [HI EJ E] , proposed by one 
of the authors of this letter: N.S.M.B.. 

The origin of families is not understood [17] so far. There are several attempts in the 
literature trying to understand the origin of families. All of them, however, in one or another 
way (for example through choices of appropriate groups) simply postulate that there are at 
least three families, as does the standard model. Proposing the (right) mechanism for gen- 
erating families is to our understanding one of the most promising guides to physics beyond 
the standard model. The approach unifying spin and charges is offering the mechanism for 
the appearance of families. It introduces the second kind [HI EJ [7J [9] of the Clifford algebra 
objects, which generates families as the equivalent representations to the Dirac spinor repre- 
sentation [IB] • The approach predicts from the simple starting action for the fermions, which 
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carry two kinds of the Clifford algebra objects' quantum numbers, more than the observed 
three families. It predicts two times four families with masses several orders of magnitude 
bellow the unification scale of the three observed charges. Since due to the approach (af- 
ter assuming a particular way of nonperturbative breaking the starting symmetry) the fifth 
family decouples in the Yukawa couplings from the lower four families [8], the quarks and 
the leptons of the fifth family are stable as required by the condition iii.. Since the masses 
of all the members of the fifth family lie much above the known three and the predicted 
fourth family masses (the fourth family might according to the first very rough estimates be 
even seen at the LHC), the baryons made out of the fifth family form small enough clusters, 
so that their scattering amplitude among themselves and with the ordinary matter is small 
enough and also the number of clusters is low enough to fulfil the conditions i. and iii.. 

There are several assessments about masses of a possible (non stable) fourth family of 
quarks and leptons, which follow from the analyses of the existing experimental data and the 
cosmological observations. Although most of physicists have doubts about the existence of 
more than three families, the analyses clearly show that neither the experimental electroweak 
data jH E] , nor the cosmological observations |4| forbid the existence of more than three (so 
far observed) families, as long as the masses of the fourth family quarks are higher than a few 
hundred GeV and the masses of the fourth family leptons higher than one hundred GeV. Our 
stable fifth family baryons, which might form the dark matter, also do not contradict the 
so far observed experimental data — as it is the measured (first family) baryon number and 
its ratio to the photon energy density, as long as the fifth family quarks are heavy enough 
(>1 TeV). (This is true for any stable heavy family.) Namely, all the measurements, which 
connect the baryon and the photon energy density, relate to the moment (s) in the history 
of the universe, when the baryons of the first family where formed {m\(? ~ k^T w 1 GeV) 
and the electrons and nuclei formed atoms (k^T « 1 eV). The chargeless (with respect to 
the colour and electromagnetic charges and not with respect to the weak charge) clusters 
of the fifth family were formed long before (at Tkb ~ E C5 (see Table [i])). They manifest 
after decoupling from the plasma (with their small number density and small cross section) 
(almost) only their gravitational interaction. 

In this paper we estimate the properties of the fifth family members (m 5 , <i 5 , z/ 5 , es) for 
which the approach unifying spin and charges predicts that they have the properties of 
the lower four families: the same family members with the same charges and interacting 
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correspondingly with the same gauge fields. We estimate the masses of the fifth family 
quarks, their behaviour in the evolution of the universe, their formation of clusters, properties 
of these clusters and the behaviour of clusters when scattering among themselves and with 
the ordinary matter (the first family baryons and leptons). 

We use a simple (the hydrogen-like) model [10J to estimate the size and the binding 
energy of the fifth family baryons, assuming that the fifth family quarks are heavy enough 
to interact mostly by exchanging one gluon. Solving the corresponding Boltzmann equations 
we estimate the behaviour of quarks and anti-quarks of the fifth family during the evolution 
of our universe, assuming that there is no excess of quarks over anti-quarks, concluding that 
quarks and anti-quarks, which succeeded to form neutral (with respect to the colour and 
electromagnetic charge) clusters, might now form the dark matter, while the rest disappeared 
at the colour phase transition at around 1 GeV and lower. We also estimate the behaviour 
of our fifth family clusters if hitting the DAMA/Nal— DAMA-LIBRA ^\ and CDMS [2] 
experiments estimating the limitations the DAMA/Nal experiments put on our fifth family 
quarks when recognizing that CDMS has not found any event yet. 

II. PROPERTIES OF CLUSTERS OF THE HEAVY FAMILY 

Let us assume that there is a heavy family of quarks and leptons as predicted by the 
approach unifying spin and charges, with masses several orders of magnitude greater than 
those of the known three families, decoupled in the Yukawa couplings from the lower mass 
families and with the charges and their couplings to the gauge fields of the known families 
(which all seems, due to our estimate predictions of the approach, reasonable assumptions). 
Families distinguish among themselves (besides in masses) in the family index (in the quan- 
tum number, which in the approach is determined by the second kind of the Clifford algebra 
objects operators [21 El E] S ab = |(7 a 7 6 — 7 fc 7 a ), anti-commuting with the Dirac 7 a 's), and 
(due to the Yukawa couplings) in their masses. 

For a heavy enough family the properties of baryons (protons p§ (u^u^d^), neutrons n$ 
(usdsds), A5 , A^ + ) made out of quarks M5 and d 5 can be estimated by using the non 
relativistic Bohr-like model with the - dependence of the potential between a pair of quarks 
y _ _ghc^c^ w here a c is in this case the colour coupling constant. Equivalently goes for 
anti-quarks. This is a meaningful approximation as long as the one gluon exchange is the 
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dominant contribution to the interaction among quarks, that is as long as excitations of a 
cluster are not influenced by the linearly rising part of the potential [TH]. The electromagnetic 



and weak interaction contributions are more than 10~ 2 times smaller. Which one of p$, n 5 , 
or maybe A5 or A5" 1 ", is a stable fifth family baryon, depends on the ratio of the bare 
masses m U5 and m^, as well as on the weak and the electromagnetic interactions among 
quarks. If is appropriately smaller than m Us so that the weak and electromagnetic 
interactions favor the neutron n 5 , then n 5 is a colour singlet electromagnetic chargeless 
stable cluster of quarks, with the weak charge —1/2. If m ds is larger (enough, due to the 
stronger electromagnetic repulsion among the two M5 than among the two rf 5 ) than m U5 , 
the proton p§ which is a colour singlet stable nucleon with the weak charge 1/2, needs the 
electron e§ or e± to form a stable electromagnetic chargeless cluster. An atom made out of 
only fifth family members might be lighter or not than 77,5, depending on the masses of the 
fifth family members. 

Neutral (with respect to the electromagnetic and colour charge) particles that constitute 
the dark matter can be n 5 , z/ 5 or charged baryons like p$, A5"" 1 ", Ajf , forming neutral atoms 
with e§ or e£, correspondingly. We treat the case that 77.5 as well as 77,5 form the major part 
of the dark matter, assuming that (and 77.5) are stable baryons (anti-baryons). Taking 
m U5 < m es also z/5 contributes to the dark matter. 

In the Bohr-like model we obtain if neglecting more than one gluon exchange contribution 

„ 1 /2 \ 2 m ar Tic 

2 V3 J 2 ' 5 § a c ^c 2 K ' 

The mass of the cluster is approximately m C5 c 2 ~ 3m qs c 2 (l — (| a c ) 2 ). We use the factor of | 
for a two quark pair potential and of | for a quark and an anti-quark pair potential. If treat- 
ing correctly the three quarks' (or anti-quarks') center of mass motion in the hydrogen-like 
model, allowing the hydrogen-like functions to adapt the width as presented in Appendix I, 
the factor —3 | (|) 2 \ in Eq. [l]is replaced by 0.66, and the mass of the cluster is accordingly 
3m (?5 c 2 (l — 0.22 a 2 ), while the average radius takes the values as presented in Table [TJ 

Assuming that the coupling constant of the colour charge a c runs with the kinetic energy 
—E C5 /3 and taking into account the number of families which contribute to the running 
coupling constant in dependence on the kinetic energy (and correspondingly on the mass of 
the fifth family quarks) we estimate the properties of a baryon as presented on Table [T] (the 
table is calculated from the hydrogen-like model presented in Appendix I), 
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TABLE I: The properties of a cluster of the fifth family quarks within the extended Bohr-like 
(hydrogen-like) model from Appendix I. m qs in TeV/c 2 is the assumed fifth family quark mass, 
a c is the coupling constant of the colour interaction at E ~ (—E C5 /3) (Eq{T]) which is the kinetic 
energy of quarks in the baryon, r C 5 is the corresponding average radius. Then a C5 = 7rr 2 5 is the 
corresponding scattering cross section. 

The binding energy is approximately ^ of the mass of the cluster (it is ~ The 
baryon n 5 [u^d^d^) is lighter than the baryon p 5 , (u qs d q5 d q5 ) if Amy = (m U5 — rrid 5 ) is 
smaller than (0.05,0.5,5,50,500,5000) GeV for the six values of the m q5 c 2 on Table |l( 
respectively. We see from Table [I] that the "nucleon-nucleon" force among the fifth family 
baryons leads to many orders of magnitude smaller cross section than in the case of the first 
family nucleons (<r Cg = irr 2 is from 10 _5 fm 2 for m q5 c 2 = 1 TeV to 10 _14 fm 2 for m q5 c 2 = 10 5 
TeV) . Accordingly is the scattering cross section between two fifth family baryons determined 
by the weak interaction as soon as the mass exceeds several GeV. 

If a cluster of the heavy (fifth family) quarks and leptons and of the ordinary (the lightest) 
family is made, then, since ordinary family dictates the radius and the excitation energies of 
a cluster, its properties are not far from the properties of the ordinary hadrons and atoms, 
except that such a cluster has the mass dictated by the heavy family members. 

III. DYNAMICS OF A HEAVY FAMILY BARYONS IN OUR GALAXY 

There are experiments [HE] which are trying to directly measure the dark matter clusters. 
Let us make a short introduction into these measurements. We shall treat our fifth family 
clusters in particular. The density of the dark matter p dm in the Milky way can be evaluated 
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from the measured rotation velocity of stars and gas in our galaxy, which appears to be 
approximately independent of the distance r from the center of our galaxy. For our Sun this 
velocity is Vs ~ (170 — 270) km/s. pd m is approximately spherically symmetric distributed 
and proportional to -\. Locally (at the position of our Sun) pd m is known within a factor 
of 10 to be po ~ 0.3 GeV/(c 2 cm 3 ), we put pd m = Po£ P , with | < e p < 3. The local velocity 
distribution of the dark matter cluster Vdmi, in the velocity class i of clusters, can only be 
estimated, results depend strongly on the model. Let us illustrate this dependence. In a 
simple model that all the clusters at any radius r from the center of our galaxy travel in all 
possible circles around the center so that the paths are spherically symmetrically distributed, 
the velocity of a cluster at the position of the Earth is equal to vs, the velocity of our Sun 
in the absolute value, but has all possible orientations perpendicular to the radius r with 
equal probability. In the model that the clusters only oscillate through the center of the 
galaxy, the velocities of the dark matter clusters at the Earth position have values from zero 
to the escape velocity, each one weighted so that all the contributions give pd m - Many other 
possibilities are presented in the references cited in [TJ. 

The velocity of the Earth around the center of the galaxy is equal to: ve = v,s + ves, 
with ves — 30 km/s and v v s s ^ ~ cos 8 smut, 9 = 60°. Then the velocity with which the 
dark matter cluster of the i- th velocity class hits the Earth is equal to: VdmEi = Vdmi — ve- 
uj determines the rotation of our Earth around the Sun. 

One finds for the flux of the dark matter clusters hitting the Earth: $ rfm = J2i 7^ \vdmi~ 
ve\ to be approximately (as long as \^ d VE t^ s \ is small) equal to 

1 \ Pdmi r 1 — * —* 1 -, ^drrii Vs t r n \ 

*rfm ~ 2^ {\Vdmi ~ v s \ - v ES ■ ZTi}- (2) 

i rn C5 \Vdmi - v s \ 

Further terms are neglected. We shall approximately take that J2i \vdm% — vs\p&mi ~ 
£v dmS £ P vsPo, and correspondingly J2i $es • ^Z^Z] ~ v ES^v dmS cos 6 sincjt, (determining 
the annual modulations observed by DAMA [1]). Here | < e VdmS < 3 and 1 < < 3 

''' v dmS 

are estimated with respect to experimental and (our) theoretical evaluations. 

Let us evaluate the cross section for our heavy dark matter baryon to elastically (the 
excited states of nuclei, which we shall treat, I and Ge, are at ~ 50 keV or higher and are 
very narrow, while the average recoil energy of Iodine is expected to be 30 keV) scatter on 
an ordinary nucleus with A nucleons a a = -|s < |^c 5 a| > 2 in\. For our heavy dark matter 
cluster is rriA approximately the mass of the ordinary nucleus. In the case of a coherent scat- 
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tering (if recognizing that A = ^ is for a nucleus large enough to make scattering coherent 
when the mass of the cluster is 1 TeV or more and its velocity ~ v$), the cross section is 
almost independent of the recoil velocity of the nucleus. For the case that the "nuclear force" 
as manifesting in the cross section 7r (r C5 ) 2 in Eq.([T]) brings the main contribution [20J the 
cross section is proportional to (3A) 2 (due to the square of the matrix element) times (A) 2 
(due to the mass of the nuclei rriA ~ 3Am qi , with m qi c 2 ~ 3 )■ When m q5 is heavier than 
10 4 TeV/c 2 (Table [I]), the weak interaction dominates and a a is proportional to (A — Z) 2 A 2 , 
since to Z° boson exchange only neutron gives an appreciable contribution. Accordingly 
we have, when the "nuclear force" dominates, a a ~ cr A A e a , with which is nr 2 5 s anucl 

and with ^ < e a , < 30. e a , takes into account the roughness with which we treat 

30 u nucl u nucl o 

our heavy baryon's properties and the scattering procedure. When the weak interaction 
dominates, e a is smaller and we have a e a = ( m ^ F ^-) 2 e aweak (= (1CT 6 im) 2 e aweak ), 
i < £a weak < 1- The weak force is pretty accurately evaluated, but the way how we are 
averaging is not. 

IV. DIRECT MEASUREMENTS OF THE FIFTH FAMILY BARYONS AS DARK 
MATTER CONSTITUENTS 

We are making very rough estimations of what the DAMA pQ and CDMS j2] experiments 
are measuring, provided that the dark matter clusters are made out of our (any) heavy family 
quarks as discussed above. We are looking for limitations these two experiments might put 
on properties of our heavy family members. We discussed about our estimations and their 
relations to the measurements with R. Bernabei [T3] and J. Filippini [13J. Both pointed 
out (R.B. in particular) that the two experiments can hardly be compared, and that our 
very approximate estimations may be right only within the orders of magnitude. We are 
completely aware of how rough our estimation is, yet we conclude that, since the number 
of measured events is proportional to (m C5 )~ 3 for masses ~ 10 4 TeV or smaller (while for 
higher masses, when the weak interaction dominates, it is proportional to (m C5 ) _1 ) that even 
such rough estimations may in the case of our heavy baryons say whether both experiments 
do at all measure our (any) heavy family clusters, if one experiment clearly sees the dark 
matter signals and the other does not (yet?) and we accordingly estimate the mass of our 
cluster. 
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Let Na be the number of nuclei of a type A in the apparatus (of either DAMA [T], which 
has 4 ■ 10 24 nuclei per kg of I, with Aj = 127, and Na, with A^ a = 23 (we shall neglect Na), 
or of CDMS |2j, which has 8.3 • 10 24 of Ge nuclei per kg, with Aq z ~ 73). At velocities of a 
dark matter cluster VdmE ~ 200 km/s are the 3A scatterers strongly bound in the nucleus, so 
that the whole nucleus with A nucleons elastically scatters on a heavy dark matter cluster. 
Then the number of events per second (Ra) taking place in Na nuclei is due to the flux 
$dm and the recognition that the cross section is at these energies almost independent of 
the velocity equal to 

Ra = N A ^ a {A) v s e VdmS e p (1 + v -™ cosfl sinu,t). (3) 

m c 5 ' £ v dmS V S 

Let ARa mean the amplitude of the annual modulation of Ra 

AR a = R A (uJt = J) - R A (oot = 0) = N A R A* — cos#, (4) 

where Ro = cr Vs e, -Ro is for the case that the "nuclear force" dominates Ro ~ 
( — 3hc ) 2 -Po_ V g £ with e = e e v , „ Q e a ,. Ro is therefore proportional to m~ 3 . We esti- 
mated 10 -3 < e < 10 3 , which demonstrates both, the uncertainties in the knowledge about 
the dark matter dynamics in our galaxy and our approximate treating of the dark matter 
properties. (When for m qs c 2 > 10 4 TeV the weak interaction determines the cross section 
Rq is in this case proportional to ^V*-) We estimate that an experiment with Na scatterers 
should measure the amplitude RA£ cu tA, with e cu tA determining the efficiency of a particular 
experiment to detect a dark matter cluster collision. For small enough £vdmES Ess cos q we 
have 

Ra £ C utA ~ N a Ro A 4 e cut A = AR A e cut A VdmS — tt- (5) 

£ VdmES V ES COS0 

If DAMA [T] is measuring our heavy family baryons scattering mostly on / (we neglect the 
same number of Na, with A = 23), then the average Rj is 

Dp ^ A D £ VdmS V S / R \ 

£ v dmES V ES cos 60 

with ARdama ~ ARi Ecutdama, this is what we read from their papers pQ. In this rough 
estimation most of unknowns about the dark matter properties, except the local velocity of 
our Sun, the cut off procedure (s cu tdama) and £vdmS (estimated to be \ < £vdmS < 3), are 

£v dmES 6 £v dmES 

hidden in AR dama . If we assume that the Sun's velocity is v$ = 100, 170, 220, 270 km/s, we 



find v^cos e = ^> ^> res P ec tively. (The recoil energy of the nucleus A = I changes cor- 
respondingly with the square of v s .) DAMA/Nal, DAMA/LIBRA [I] publishes AR dama = 

0.052 counts per day and per kg of Nal. Correspondingly is Ri e cut dama = 0, 052 £vdmS 

counts per day and per kg. CDMS should then in 121 days with 1 kg of Ge (A = 73) detect 
RGeBcutcdms ~ M (^l) 4 £ ™ tcdms I^s Q Q52 . m eve nts, which is for the above 

Q.yJ 1Z / ^cut dama £» Jm p, vse cos V 

measured velocities equal to (10, 16,21,25) £e " tcdms £vdmS . CDMS M has found no event. 

The approximations we made might cause that the expected numbers (10,16,21,25) 
multiplied by 6cutGe £ " dmS are too high (or too low!!) for a factor let us say 4 or 10. If in 

Scut I £ v dmES 

the near future CDMS (or some other experiment) will measure the above predicted events, 
then there might be heavy family clusters which form the dark matter. In this case the 
DAMA experiment puts the limit on our heavy family masses (Eq.(j5])). 

Taking into account all the uncertainties mentioned above, with the uncertainty with 
the "nuclear force" cross section included (we evaluate these uncertainties to be 10~ 4 < 
e < 3 • 10 3 ), we can estimate the mass range of the fifth family quarks from the DAMA 
experiments: (m q5 c 2 ) 3 = ^^Nj A 4 tt (^f) 2 p c 2 v ES cosOe = (0.3 • 10 7 ) 3 e (g) 2 GeV. 
The lower mass limit, which follows from the DAMA experiment, is accordingly m q5 c 2 > 200 
TeV. Observing that for m q5 c 2 > 10 4 TeV the weak force starts to dominate, we estimate 
the upper limit m q5 c 2 < 10 5 TeV. Then 200 TeV < m q& c 2 < 10 5 TeV. 



V. EVOLUTION OF THE ABUNDANCE OF THE FIFTH FAMILY MEMBERS 
IN THE UNIVERSE 

To estimate the behaviour of our stable heavy family quarks and anti-quarks in the 
expanding universe we need to know: 

i. ) the masses of our fifth family members, 

ii. ) their particle — anti-particle asymmetry, 

iii. ) their thermally averaged scattering cross sections (as the function of the temperature) for 
scattering iii. a.) into all the relativistic quarks and anti-quarks of lower families (< av > q q), 
iii.b.) into gluons (< av > gg ), iii.c.) into (annihilating) bound states of a fifth family quark 
and an anti-quark (< av >( q q) b ), iii.d.) into bound states of two fifth family quarks and 
into the fifth family baryons (< av > C5 ) (and equivalently into two anti-quarks and into 
anti-baryons), 
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iv.) the probability for quarks and anti-quarks of the fifth family to annihilate at the colour 
phase transition (T/c fe « 1 GeV). 

The quarks and anti-quarks start to freeze out when the temperature of the plasma falls 
close to m qs c 2 /kb is the Boltzmann constant). They are forming clusters (bound states) 
when the temperature falls close to the binding energy. When the three quarks or three 
anti-quarks of the fifth family form a colourless baryon (or anti-baryon) , they decouple from 
the rest of the plasma due to small scattering cross section manifested by the average radius 
presented in Table [1} 

We assume in this paper that there is no asymmetry between quarks and anti-quarks of 
the fifth family. One evaluates that at the colour phase transition (Tkb ~ lGeV) the ratio 
of the scattering time between two coloured quarks (of any kind) and the Hubble time is 
of the order of ~ 10~ 18 . Accordingly, although the number of the fifth family quarks and 
anti-quarks is of the order of 10~ 13 smaller than the number of the quarks and anti-quarks 
of the rest of families (as show the solutions of the Boltzmann equations presented bellow), 
the fifth family quarks and anti-quarks have enough opportunity during the expansion time 
(from 10~ 7 s to a few seconds) to deplete completely. The same would happen to all the 
lower families' quarks and anti-quarks (going due to the Yukawa couplings to the first family 
members), if there would be no quark — anti-quark asymmetry. 

To see how many fifth family quarks and anti-quarks succeed to form the fifth family 
baryons and anti-baryons we must solve the Boltzmann equations as a function of time (or 
temperature). Since we do not know the mass of the fifth family members (the estimations 
from the approach unifying spin and charges predict that it must be higher than a few TeV 
and lower than, say, a few 10 5 TeV), we take it as a parameter. The interaction due to 
one gluon exchange is dominant among so massive fifth family members and it is at the 
same time also much larger than the weak and the electromagnetic interaction. Since the 
one gluon exchange is (up to the group properties and the coupling constants) equivalent 
to the one photon exchange, we use for the cross sections (cited above) the equivalent cross 
sections from the electromagnetic case. The fifth family quark mass follows from comparing 
the calculated fifth family baryon and anti-baryon number density multiplied by the mass 
of the clusters with the today's dark matter density. 

We follow (as much as possible), when estimating the number density of the fifth family 
quarks n qs and anti-quarks rig 5 clustered into baryons (with the number density n C5 ) and 
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anti-baryons (ne 5 ), which to our prediction form the dark matter today, the ref. [3], chapter 
3. n qs is the number density of all the fifth family quarks of any colour and spin and 
correspondingly is assumed for the other number densities. The following cross sections are 
needed in the Boltzmann equations 



where v is the relative velocity between the fifth family quark and its anti-quark, or between 



averaged scattering cross section times the relative velocity: i. < ov > qq - for all the pairs of 
the fifth family quarks and anti-quarks into all the lower mass (of the four families') quarks 
and anti-quarks, which are, while scattering takes place, ultra relativistic. ii. < ov > gg 
for scattering into gluons. iii. < ov > C5 for two quarks (or two anti-quarks) to scatter 
into a bound state of two quarks (anti-quarks) and from two to three quarks (anti-quarks) 
colourless clusters. We use the equivalent expression as for scattering of an electron and a 
proton into the bound state of a hydrogen. The parameter t] cs takes care of scattering of 
two quarks (anti-quarks) into three colourless quarks (or anti-quarks), which are the fifth 
family baryons (anti-baryons). iv. < ov >( qq ) b for scattering into a bound state of the fifth 
family quark and anti-quark, which annihilate in the time T( q q) b < 10~ 28 s. t]( q q) b takes care 
of the roughness of the used formula. 

ot is the Thompson-like scattering cross section of gluons on quarks (or anti-quarks), 
responsible for destroying the bound states of baryons. 

Let T be the today's black body radiation temperature, T(t) the actual (studied) temper- 
ature, a 2 (T°) = 1 and a 2 (T) = a 2 (T(t)) is the metric tensor component in the expanding flat 
universe — the Friedman-Robertson- Walker metric: diagg^ = (1, — a(t) 2 , — a(t) 2 , — a(t) 2 ), 
(f) 2 = ^Y~P-> w ith p = f| g* T 4 , T = T(t), g* measures the number of degrees of freedom of 




(7) 



two quarks, E C5 is the binding energy for a cluster (Eq. [I]). < ov > is the thermally 
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those of the four family members (f) and gauge bosons (b), which are at the treated tem- 
perature T ultra-relativistic (g* = J2ieh 9% + | Sief 9i)- Ho ~ 1-5 1CT 42 is the present 
Hubble constant and G = t^t, m^c 2 rs 1.2 • 10 19 GeV. 

Let us write down the Boltzmann equation, which treats in the expanding universe the 
number density of all the fifth family quarks as a function of time t. The fifth family quarks 
scatter with anti-quark into all the other relativistic quarks and anti-quarks (< ov > q q) and 
into gluons (< ov > gg )- At the beginning, when the quarks are becoming non-relativistic and 
start to freeze out, the formation of bound states is negligible. One finds [3] the Boltzmann 
equation for the fifth family quarks n q5 (and equivalently for anti-quarks n q5 ) 

3 d(q% 5 ) ( )_(o) ( n qT n qb n q n- q 

dt - <GV >q * "* ^ { r4?4? nS°> 4 0) 

<™> 99 <^(-lf% + J ^)- (8) 
\ n\ b 'n\ 5 > nynyj 

Let us tell that nf ] = # ( "gffi* )^"^" for m iC 2 » Tk b and to ^ (^) 3 for m iC 2 « 
Tkb. Since the ultra-relativistic quarks and anti-quarks of the lower families are in the 
thermal equilibrium with the plasma and so are gluons, it follows mwo) = 1 = wTm ■ 
Taking into account that (aT) 3 g*(T) is a constant it is appropriate [3] to introduce a new 

2 

parameter x = m £ 6 £ and the quantity Y q5 = n qs (^) 3 , = n q °J (^§0 3 - When taking 
into account that the number of quarks is the same as the number of anti-quarks, and that 



- - with h m = J^f Eq. 8 transforms into ^ = ^ (Vf ) 2 - Y 2 ), with 



dt ' x ' m V 45 hcrripic 2 '""II ^ x x 95 

A gs = ( < ^ > gg+ < ^>gg) m g5 c _ j g fais equation which we are solving (up to the region of x 
when the clusters of quarks and anti-quarks start to be formed) to see the behaviour of the 
fifth family quarks clS cL function of the temperature. 

When the temperature of the expanding universe falls close enough to the binding energy 
of the cluster of the fifth family quarks (and anti-quarks), the bound states of quarks (and 
anti-quarks) and the clusters of fifth family baryons (in our case neutrons n§) (and anti- 
baryons 77.5 — anti-neutrons) start to form. The corresponding Boltzmann equation for the 
number of baryons n C5 reads 

dt <<JV> <* "* [{ n (o)) n (o)j- W 

Introducing again Y C5 = n C5 (^) 3 , = ng> (^) 3 and A C5 = " ' with the same 
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x and h m as above, we obtain the equation ^j 5 - = ^ (Y 2 — Y" C5 Y 



Co) 



2 CK 



The number density of the fifth family quarks n q5 (and correspondingly Y q5 ), which has 
above the temperature of the binding energy of the clusters of the fifth family quarks (almost) 
reached the decoupled value, starts to decrease again due to the formation of the clusters of 
the fifth family quarks (and anti-quarks) as well as due to forming the bound state of the 
fifth family quark and an anti-quark, which annihilates into gluons. It follows 



dt ^»»^c 5 " qs >*q B I (0) I ^ (0) n I (0) 



+ 



,(o)Jo) / _ V%. , n q n q | 

-qq '^5 "- 95 (0) (0) ^ (0) (0) 1 

' n q5 'n y q 5 ' n q 'n\ 

M (0) / V% , ra g"g j , , ( V , 

'99 "-g 5 ""95 (0) (0) ^ (0) (0) I ' - l ■ ' 

1 n q5 'nq 5 ' n g 'rig 



with and ?? C5 defined in Eq. [7j Introducing the above defined Y qs and F C5 the Eq. 10 
transforms into ^ = ^ (-V, 2 +F5 ^£) + ^ (-Y*) + ^ «°> 2 -Y ? 2 ), * A (9j06 = 



<( ™k (g (^ e )3 g5C (and with the same x and h m as well as A C5 and X q5 as defined above). We 
solve this equation together with the above equation for Y C5 . 

Let us look at also at the Thompson scattering of gluons on the bound states, destroying 
clusters, which starts to be negligible when the rate for gluons to scatter off the quarks 



(n q5 a T c) starts to be smaller than the expansion rate (H = \f Sn ^* ^^J ), with g* defined 
above. Recognizing that the binding energy of Table [i] is approximately ^ the mass of 
the fifth family quarks we get the requirement that the bound states can be formed when 



n. 



1 

,, 5 << 3.10~ 25 ( ^y ) 4 fm~ 3 , which for m q5 c 2 = 1 TeV gives n q5 « 3.10" 13 fm~ 3 and for 
m qs c 2 = 10 TeV gives n q5 « 3.10 -9 fm~ 3 . One can easily check from the solutions of the 
Boltzmann equations that this requirements are fulfilled. 

Solving the two Boltzmann equations (Eqs. |8j [9]) we obtain the number density of the 
fifth family quarks n qs (and anti-quarks) and the number density of the fifth family baryons 
n C5 (and anti-baryons) as a function of the mass of the fifth family quarks and the two 
parameters t] C5 and f]( qq ) 3 - The evaluations are made, as we explained above, within the 
approximation that the formation of the baryons is determined by the formation of the 
two quark (anti-quarks) clusters, corrected by the parameter r] C5 (Eq. [7]). The parameter 
rjr qq -\ b corrects the formation of a quark — anti-quark pair and its annihilation into gluons. In 
Diagram V both number densities (multiplied my (^) 3 , which is Y q5 and Y C5 , respectively 
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for the quarks and the clusters of quarks) as a function of m ^ for ?7( g g) 3 = 1 and r] C5 = ^ 
are presented. The calculation is performed up to Tk b = 1 GeV. 




FIG. 1: The dependence of the two number densities n q5 (of the fifth family quarks) and n C5 (of the 

2 

fifth family clusters) as the function of is presented for the special values m q5 c 2 = 71 TeV, 

Vcr, = and V(qQ)b = L We take 9* = 91-5. 

2 

The quarks and anti-quarks are at high temperature ( "^ << 1) in thermal equilibrium 
with the plasma (as are also all the other families and bosons of lower masses). As the 
temperature of the plasma (due to the expansion of the universe) drops close to the mass of 
the fifth family quarks, quarks and anti-quarks scatter into all the other (ultra) relativistic 
fermions and bosons. At the temperature close to the binding energy of the quarks in 
a cluster, the clusters of baryons start to be formed. We evaluated the number density 
n qs (T) (^) 3 = Y qs of the fifth family quarks (and anti-quarks) and the number density of 
the fifth family baryons n Cs (T) (^) 3 = Y C5 for several choices of m q6 ,r) cs and Tj^ b up to 
Tk b = 1 GeV = 

/,From the calculated decoupled number density of baryons and anti-baryons of the fifth 
family quarks (and anti-quarks) n C6 (Ti) at temperature T\kb = 1 GeV, where we stopped 
our calculations as a function of the quark mass and of the two parameters r) CB and fj( qq ) b , 
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m qr <? 
TeV 


Viqqh = 1 


T7 — 1 

'(99)6 3 


^(^t. = 3 


V(qq) b = 10 


Vc 5 = 1 


19 


11 


37 




Vc 5 = 3 


15 


9.5 


27 




Vcs = I 


25 


14 


54 




Vet = 10 


13 




22 




^5 = 10 


39 


20 


84 




^cs = 50 


71 






417 



TABLE II: The fifth family quark mass is presented, calculated for different choices of rj C5 (which 
takes care of the probability that a colourless cluster of three quarks (anti-quarks) instead of two are 
formed) and of 'H(qq) b (which takes care of the annihilation of a bound state of quark — anti-quark) 



from Eqs. (|ll| |9j |8j) . 

the today's mass density of the dark matter follows 

Pdm = QdmPcr = 2 m C5 7l Cs {T\ 



(11) 



TJ g*(T y 

with T = 2.5 -10- 4 f , g*(T ) =2 + | -3 ■ (^) 4 / 3 , g*{T x ) = 2 + 2 ■ 8 + |(5 • 3 ■ 2 ■ 2 + 6 ■ 2 • 2) 
and p cr c 2 ~ 3 87r °g ~ 5.7 ■ lO 3 ^ , factor 2 counts baryons and anti-baryons (since the spin 
of baryons is taken into account in n CB ). 

The influence of the choice of the parameters 7] C5 and t]( q q) b on the number of the fifth 
family baryons and anti-baryons (together with the evaluation of the numerical errors) is 
used as a measure for the accuracy with which we evaluated the fifth family mass. 

We read from Table [IT] the mass interval for the fifth family quarks' mass, which fits 



Eqs. (11 9 8: 



10 TeV < m q5 c 2 < a few ■ 10 2 TeV. 



(12) 



From this mass interval we estimate from Table |T] the cross section for the fifth family 
neutrons 7r(r C5 ) 2 : 



l(T 8 fm 2 < a C5 < 10~ 6 fm 2 . 
(It is at least 10~ 6 smaller than the cross section for the first family neutrons. 



(13) 
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VI. CONCLUDING REMARKS 



We estimated in this paper the possibility that a new stable family, predicted by the 
approach unifying spin and charges [13 El IB] , forms baryons which are the dark matter 
constituents. The approach (proposed by S.N.M.B.) is to our knowledge the only proposal 
in the literature so far which offers the mechanism for generating families, if we do not count 
those which on one or another way just assume more than three families. We evaluated the 
limits on the properties of the stable fifth family quarks due to the cosmological observations 
and the direct experiments. 

We use the simple hydrogen-like model to evaluate the properties of these heavy baryons 
and their interaction among themselves and with the ordinary nuclei. We take into account 
that for masses of the order of 1 TeV/c 2 or larger the one gluon exchange determines the force 
among the constituents of the fifth family baryons. Studying the interaction of these baryons 
with the ordinary matter we take into account that the weak interaction starts to dominate 
over the "nuclear interaction" which the fifth family neutron manifests for massive enough 
clusters (m qs > 10 4 TeV), while the non relativistic fifth family baryons interact among 
themselves with the weak force only. 

We assume that in the evolution of our universe q$ and were formed with no asymme- 
try. We study the freeze out procedure of the fifth family quarks and anti-quarks and the 
formation of baryons and anti-baryons up to the temperature Tkf, = 1 GeV, when the colour 
phase transition starts which depletes all the fifth family quarks and anti-quarks while the 
colourless fifth family neutrons with very small scattering cross section decouples long before 
(at Tk b = 100 GeV). 

While the measured density of the dark matter does not put much limitation on the 
properties of heavy enough clusters, the DAMA experiments [1] limit (provided that they 
measure our heavy fifth family clusters) the quark mass to: 200 TeV < m q5 c 2 < 10 5 TeV. 
The estimated cross section for the dark matter cluster to (elastically, coherently and non- 
relativisically) scatter on the (first family) nucleus is in this case determined on the lower 
mass limit by the "nuclear force" of the fifth family clusters ((3 ■ 10~ 5 A 2 fm) 2 ) and on the 
higher mass limit by the weak force ((A(A — Z) 10 _6 fm) 2 ). 

The cosmological evolution suggests for the mass limits the range 10 TeV < m q5 c 2 < 
a few • 10 2 TeV and for the scattering cross sections 10~ 8 fm 2 < a C5 < 10~ 6 fm 2 . Accordingly 
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we conclude that if the DAMA experiments are measuring our fifth family neutrons, the 
mass of the fifth family quarks is a few hundred TeV /c 2 . 

In the ref. [TB] [2TJ the authors study the limits on a scattering cross section of a heavy 
dark matter cluster with the ordinary matter. They assume (approximately) the same 
number of particles and antiparticles in the dark matter. They treat the conditions under 
which would the heat flow which would follow from the annihilation of dark matter particles 
and anti-particles in the Earth core start to be noticable. Using their limits we conclude 
that our fifth family baryons of the mass of a few hundreds TeV/c 2 have for a factor more 
than 100 too small scattering amplitude with the ordinary matter to cause a measurable 
heat flux. 

Our rough estimations predict that, if the DAMA experiments observe the events due to 
our (any) heavy family members, (or any heavy enough family cluster with small enough 
cross section), the CDMS experiments [2] will in the near future observe a few events as 
well. If CDMS will not confirm the heavy family events, then we must conclude, trusting 
the DAMA experiments, that either our fifth family clusters have much higher cross section 
due to the possibility that is lighter than so that their velocity slows down when 
scattering on nuclei of the earth above the measuring apparatus bellow the threshold of the 
CDMS experiment (and that there must be in this case the fifth family quarks — anti-quarks 
asymmetry) while the DAMA experiment still observes them, or the fifth family clusters 
(any heavy stable family clusters) are not what forms the dark matter. 

Let us comment again the question whether there is at all possible (due to electroweak 
experimental data) that there exist more than three up to now observed families, that is, 
whether the approach unifying spin and charges by predicting the fourth and the stable fifth 
family (with neutrinos included) contradict the observations. In the ref. [15] the properties 
of all the members of the fourth family were studied (for one particular choice of breaking the 
starting symmetry). The predicted fourth family neutrino mass is at around 100 GeV/c 2 or 
higher, therefore it does not due to the detailed analyses of the electroweak data done by the 
Russian group [H] contradict any experimental data. (The stable fifth family neutrino has 
due to our calculations considerably higher mass. Accordingly none of these two neutrinos 
contradict the electroweak data. They also do not contradict the nucleosynthesis, since to 
the nucleosynthesis only the neutrinos with masses bellow the electron mass contribute. 
The fact that the fifth family baryons might form the dark matter does not contradict the 
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measured (first family) baryon number and its ratio to the photon energy density as well, 
as long as the fifth family quarks are heavy enough (>1 TeV). All the measurements, which 
connect the baryon and the photon energy density, relate to the moment (s) in the history of 
the universe, when the baryons (of the first family) where formed (rriic 2 ~ k^T = 1 GeV and 
lower) and the electrons and nuclei were forming atoms (fc&T w 1 eV). The chargeless (with 
respect to the colour and electromagnetic charges, not with respect to the weak charge) 
clusters of the fifth family were formed long before (at Tkf, « E C5 (Table [I])). They manifest 
after decoupling from the plasma (with their small number density and small cross section) 
(almost) only their gravitational interaction. 

Let us conclude this paper with the recognition: If the approach unifying spin and charges 
is the right way beyond the standard model of the electroweak and colour interaction, then 
more than three families of quarks and leptons do exist, and the stable (with respect to the 
age of the universe) fifth family of quarks and leptons is the candidate to form the dark 
matter. 



VII. APPENDIX I. THREE FIFTH FAMILY QUARKS' BOUND STATES 

We look for the ground state solution of the Hamilton equation H\ip) = E C5 for a 
cluster of three heavy quarks with 

3 r> 2 2 3 he n 

H = Y.^zr-{ E uF^> ( 14 ) 



i=i Zrriqs 3 j<3 - =1 |ij Xj ] 



in the center of mass motion 



- - -, -, -> Xl+%2 £ X-i + X 2 + £3 

x = x 2 -X!, y = x 3 — , R= , (15) 

assuming the antisymmetric colour part {\^) c ,a)i symmetric spin and weak charge 
part GV^iuspin^) and symmetric space part (\if)) spacej 5). For the space part we 
take the hydrogen-like wave functions ipa(x) =^7=f e~' x ' //a and ipb(y) = ^/=p e~^^ b , 
allowing a and b to adapt variationally. Accordingly (xx, x*2, a^lV^spaceS = 
Af ( , ip a (x)'ipb(y) + symmetric permutations). It follows (^i,^2,^3|'0)space5 = 

M (2^ a (x)Mv) + 2Mv- l)Ml + f )) + ^a{y + DMl 



i 



,,2 



The Hamiltonian in the center of mass motion reads H = , - H ^ h 

2 (_|5.) 2(=^L) 2 - 3m -J5 

I Tic a c 1 1 + + ■ Varying the expectation value of the Hamiltonian with respect 
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to a and b it follows: f = 1.03, aa ^ c = 1.6. 

Accordingly we get for the binding energy E Cs = 0.66 m qs c 2 a 2 and for the size of the 
cluster J(\x 2 — Xi| 2 ) = 2.5 — — — . 

To estimate the mass difference between w 5 and d§ for which u^,d 5 d 5 is stable we treat 
the electromagnetic (a e im) and weak (q; w ) interaction as a small correction to the above 
calculated binding energy: H' = a e i mw he ( - + r— ^ + ). a e i mw stays for electro- 
magnetic and weak coupling constants. For m q5 = 200 TeV we take a e i mw = t^j, then 
|m U5 - m ds | < § ^ C5 (fae g m)2 = 0.5 ■ 10- 4 m q5 c 2 . 
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